arXiv:1504.01256vl [nlin.AO] 6 Apr 2015 


Non-linear Dynamics, Emergent Behaviors and Controlled Expansions: Towards 

Effective Modeling of the Congested Traffic 


Bo Yang^, Xihua Xu^’^, John Z.F. Pang^, and Christopher Monterola^ 

^ Complex Systems Group, Institute of High Performance Computing, A*STAR, Singapore, 138632. 

^ Department of Mathematics, National University of Singapore, 119076, Singapore, and 
^ Beijing Computational Science Research Center, Beijing 100084, PR China. 

(Dated: April 7, 2015) 

We propose a framework for constructing microscopic traffic models from microscopic acceleration 
patterns that can in principle be experimental measured and proper averaged. The exact model 
thus obtained can be used to justify the consistency of various popular models in the literature. 
Assuming analyticity of the exact model, we suggest that a controlled expansion around the constant 
velocity, uniform headway “ground state” is the proper way of constructing various different effective 
models. Assuming a unique ground state for any fixed average density, we discuss the universal 
properties of the resulting effective model, focusing on the emergent quantities of the coupled non¬ 
linear ODEs. These include the maximum and minimum headway that give the coexistence curve 
in the phase diagram, as well as an emergent intrinsic scale that characterizes the strength of 
interaction between clusters, leading to non-trivial cluster statistics when the unstable ground state 
is randomly perturbed. Utilizing the universal properties of the emergent quantities, a simple 
algorithm for constructing an effective traffic model is also presented. The algorithm tunes the 
model with statistically well-defined quantities extracted from the flow-density plot, and the resulting 
effective model naturally captures and predicts many quantitative and qualitative empirical features 
of the highway traffic, especially in the presence of an on-ramp bottleneck. The simplicity of the 
effective model provides strong evidence that stochasticity, diversity of vehicle types and modeling of 
complicated individual driving behaviors are not fundamental to many observations of the complex 
spatiotemporal patterns in the real traffic dynamics. We also propose the nature of the congested 
phase can be well characterized by the long lasting transient states of the effective model, from 
which the wide moving jams evolve. 

PACS numbers: 89.40.+k, 47.54.-|-r, 64.60.Cn, 64.60.Lx 


I. INTRODUCTION 

Modeling the dynamics of the highway traffic flow has 
been the endeavor of researchers in many disciplines for 
the last fifty years [iHl]. Various different models have 
been proposed to describe both the free and the con¬ 
gested phase of the traffic flow[5Ml irUHTH ITMT5] . Most 
of these models can capture the low density free flow 
phase, and the wide moving jams when the density is 
high. Kerner[2n] first suggested that some essential em¬ 
pirical features are not captured by most of these models; 
there exists a “synchronized phase” that can be distin¬ 
guished by a scattering of data points covering a two- 
dimensional region on the flow-density plane. It is argued 
this phase is qualitatively different from the wide mov¬ 
ing jams, particularly when a bottleneck at the highway 
is present pri 1^ . This raises the questions of the rele¬ 
vance of the popular general motor (GM) model classes 
to real traffic systems, because these models only describe 
a two-phase transition|20j. 

Most of the three-phase models are constructed by 
putting in a “synchronization gap” by hand[T(71 ITTl 1^ . 
at the cost of making the models more sophisticated 
with significantly more parameters. These, together with 
other three-phase models [TS], reproduce the “synchro¬ 
nized phase” with a multitude of steady states in the 
congested phase. However, Helbing et.al[27] pointed out 
that the characterization of the complex congested states 


of the traffic flow as a single synchronized phase is del¬ 
icate, and with properly adjusted parameters some GM 
models can reproduce many empirical observations at the 
highway bottleneck[53. One should also note inhomoge¬ 
neous road conditions and vehicle types [28], as well as 
stochastic driving behaviors can contribute to the scat¬ 
tering of the flow-density plot. More importantly, it is 
not well understood if the empirical data in the congested 
phase comes from equilibrium/steady traffic conditions, 
or from slowly evolving transient ones. 

Unlike conventional physical systems, where the con¬ 
struction of the physical models are guided and con¬ 
strained by symmetries, in traffic systems such symme¬ 
tries are conspicuously lacking. Even individual com¬ 
ponents are not identical to each other: each of them 
responds to the interaction and environment differently, 
and in a time-dependent way. While constructions of the 
traffic models are generally guided by simplicity and the 
use of physically transparent parameters, there is a cer¬ 
tain level of arbitrariness in how the model should look 
like due to the lack of a more fundamental guiding prin¬ 
ciple. This is one of the main reasons for a plethora of 
different traffic models and the controversies in the field. 

In this paper, we aim to explore such guiding princi¬ 
ples by proposing a systematic way of constructing traffic 
models, in order to remove the arbitrariness in the pos¬ 
sible forms of the models. The first part of the paper is 
to develop a general framework to obtain the master de- 
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terministic microscopic model, which can in principle be 
obtained from empirical measurements. Various simpli¬ 
fied effective models can thus be constructed via the con¬ 
trolled expansions of the master microscopic model. For 
the next part we use the published empirical data of the 
A-5 North German Highway as an example to illustrate 
how a minimal effective model should be constructed and 
tuned. We first establish a list of empirical features of the 
traffic dynamics that are commonly reported and well- 
defined. This includes most of the quantitative charac¬ 
teristics and the qualitative spatiotemporal traffic pat¬ 
terns on which the “three phase” traffic theory is based. 
This is followed by our proposal of a simple algorithm 
to construct a minimal effective model that can capture 
all these empirical features, with the intention that such 
model is the simplest possible, and it has enough predic¬ 
tive power for it to be useful for traffic optimization and 
transport engineering. 

The resulting effective model is within the frame¬ 
work of our guiding principle and is able to capture 
unambiguously the empirical features that are previ¬ 
ous thought to require the sophisticated “three phase 
models” [ini m na, or the intelligent driver models 
(IDM)[2a that have an artificial divergence when the 
headway goes to zero. More importantly our effective 
model is surprisingly simple, providing strong evidence 
that the non-linear dynamics alone is fundamental to the 
empirical features listed in Sec. A careful interpreta¬ 
tion of the emergent quantities in the system of non-linear 
ODEs is an important component of our algorithm, in 
which the parameters and functional forms in the model 
are no longer estimated, but well-defined by macroscopic 
empirical features. 


The requirement that empirical features are common 
and well-defined does not come without a price: given 
the complex nature of the traffic system, there are ob¬ 
servations of various spatiotemporal traffic patterns that 
are difficult to differentiate in a well-defined way. Thus 
the predictive power of our model, as well as our the¬ 
oretical understanding, is at best ambiguous for those 
observations. Nevertheless, we aim to achieve a clear un¬ 
derstanding of what can be unambiguously captured by 
the model and theory given the set of necessary assump¬ 
tions, and to establish a useful “reference point”: when 
better understandings of the empirical features emerge, 
one can systematically extend and generalize the model. 


The paper will be organized as follows: In Sec. we 
present a list of empirical features, defining them en¬ 
tirely from the experimental point of view, minimizing 
the tendency of pre-mature theoretical interpretations. 
In Sec. |III| we present a general framework of systemat¬ 
ically modeling the traffic dynamics from the most gen¬ 
eral equations of motion. All assumptions of the model 
are listed and connections to various microscopic mod¬ 
els in the literature are discussed. In Sec. W we study 
some of the universal properties of a class of simplified 
effective models, which belongs to the OV model class, 
and focus on the quantities emerging from the non-linear 


anisotropic interactions between vehicles. In particu¬ 
lar we introduce an emergent length scale that quanti¬ 
fies the strength of interaction between the quasisolition 
structures appearing in the solutions of the OV mod¬ 
els, leading to non-trivial cluster statistics when the con¬ 
stant headway solutions (or the ground states) of the 
OV model are randomly perturbed. This emergent scale 
plays an important role in tuning the model and defin¬ 
ing the time scale in the traffic dynamics. In Sec. 
we present a simple algorithm in constructing a mini¬ 
mal effective model based on the macroscopic empiri¬ 
cal data, and in particular their relationship with the 


emergent quantities in the traffic model, and in Sec. VI 


we present the model’s predictions of the quantitative 
empirical features, together with the characteristic spa¬ 
tiotemporal patterns, and compare them with the empir¬ 
ical observations. In Sec. IVIII we summarize our results 
and discuss the outlooks of the traffic modeling. 


II. EMPIRICAL FEATURES OF TRAFFIC 
DYNAMICS 


Many work has been published on the empirical studies 
of the highway traffic dynamics, and readers can refer to 
[11271 US] for detailed information. The reconstruction of 
the complex spatiotemporal patterns along the highway 
based on the measurement of the flow and average ve¬ 
locity at a specific location is also a subtle issue, and an 
excellent discussion can be found in|27]. One of the most 
popular techniques in analyzing the traffic data is the 
construction of the flow-density diagram, where the flow 
of traffic through a fixed cross-section is plotted against 
the average density of the vehicles on the highway. A 
schematic drawing of such plot is shown in Fig.Q. 

It is well understood that in the limit of the road den¬ 
sity p —> 0, the flow-density relationship is linear with 
very small scattering of the data. The gradient in that 
limit limp_>o dQ/dp is the maximum velocity Vmax of the 
vehicles on the highway. This is the characteristic veloc¬ 
ity of the vehicle when it is travelling freely with no other 
vehicle in the front, and is in general constrained by the 
speed limit, road conditions and the physical capabilities 
of the vehicle. 

When the density increases from zero, initially the av¬ 
erage velocity of the vehicles does not change much, as 
the vehicles hardly interact with each other. Thus the 
flow increases linearly. At intermediate densities, the 
interaction between vehicles strongly reduce the aver¬ 
age velocity, leading to a sub-linear increase of the flow 
against the density. This continues until the density 
reaches the critical density pc, at which the flow reaches a 
maximum Qmax- The maximum flow Qmax is also defined 
as the capacity of the highway. 

In general the relationship between the flow and den¬ 
sity becomes more scattered as the density increases to¬ 
wards Pc- What is more interesting is a characteristic 
discontinuous drop in Q when p increases past pc, to- 
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FIG. 1. The schematic drawing of the flow-density plot of the 
highway traffic, taken from Ref. [T]. The free flow phase is 
labeled as F, and the scattering of the flow-density plot in 
the congested phase is shaded and labeled by S. The Jam 
line labeled by J describes the wide moving jam. The density 
is normalized by the density within a wide moving jam, and 
the maximum flow Qmax is achieved at the critical density p^. 
Qdj and pdj are the flow and density of the traffic downstream 
of a wide moving jam, while the gradient of the J line gives 
the characteristic velocity of the downstream front of the wide 
moving jam. 

gether with an onset of wide scattering of the flow-density 
relationship, covering a two-dimensonal area on the flow- 
density plane. There Is clearly a phase transition here 
physically corresponding to the breakdown of the free 
flow traffic to the congested traffic. In Kerner’s “three- 
phase traffic theory” [1] , the two-dimensional scattering 
of the flow-density data points corresponds to the “syn¬ 
chronized flow”, or the “synchronized phase”. The fun¬ 
damental assumption of the “three-phase theory” is that 
each scattered data point on the flow-density plane corre¬ 
sponds to a steady state, that can be either stable or un¬ 
stable against the formation of the wide moving jam. Mi¬ 
croscopically it is postulated that the speed adaptation of 
the drivers at high vehicle density leads to a non-unique 
relationship between the average velocity and the aver¬ 
age headway. It is also postulated the two-dimensional 
covering of the “synchronized phase” is bounded by the 
free flow, and a characteristic upper and lower boundary. 

In this paper, we take a step back and stop short of 
making assumptions about the nature of the scattered 
flow-density relationship. The real highway traffic Is an 
open system and It is very difficult to verify experimen¬ 
tally if these states are steady states, or just transient 
states that last for long enough time for It to be captured 
by the sensors. Instead we just characterize the congested 
phase of the traffic dynamics by a depression of the traf¬ 
fic flow and a much larger scattering of the flow-density 
relationship as compared to the free flow when p < pc- 

There is also a well-known “hysteresis effect” m, 
whereby the traffic flow breaks down at the maximum 
flow Qmax, and after staying congested for an extended 
period of time returns back to the free flow but at a 
smaller flow. While due to the wildly fluctuating na¬ 


ture of the congested phase one cannot characterize the 
“hysteresis effect” quantitatively, one can understand it 
from the crucial empirical observation that the traffic 
flow downstream of the congested phase, and in particu¬ 
lar of the wide moving jam, is significantly smaller than 
the maximum flow Qmax[l], and this should also be pre¬ 
dicted by a useful model. 

Another important feature in the flow-density diagram 
is the “J line”, corresponding to the emergence of the 
wide moving jam in the traffic. The structure of the 
wide moving jam is surprisingly robust [1], leading to sev¬ 
eral readily measurable quantities that are characteristic 
of the traffic system. The gradient of the “J line” is ob¬ 
tained from the characteristic velocity Vj of the jam mov¬ 
ing upstream; the intersection of the “J line” with the 
free flow branch of the flow-density plot comes from the 
characteristic flow Qdj , and the density pdj of the vehicle 
right downstream of the wide moving jam. The intersec¬ 
tion of the “J line” with the density axis comes from the 
density of the vehicle, pj, within the wide moving jam, 
where all vehicles come to a stop with zero velocity. One 
should also note the following relationship 


It was noted by Kerner these characteristic quantities 
does not change with the flow and density of the vehi¬ 
cles away from the wide moving jam (both upstream and 
downstream), thus a useful traffic model should be able 
to capture these Important quantities. 

Apart from the flow density diagram, It Is also very im¬ 
portant to characterize the rich spatiotemporal patterns 
of the traffic in the congested phase. It is well established 
that the congested phase Is macroscoplcally rather homo¬ 
geneous as compared to the wide moving jams, though 
microscopically one would observe spatial and temporal 
fluctuations of both density and velocity, sometimes in 
the form of numerous narrow jams in which the vehi¬ 
cles’ velocities drop to zero momentarllv[il [^l^[?7ll^ . 
The distinction between the general congested traffic (the 
“synchronized phase” or the “general patterns”) and the 
wide moving jams Is rather difficult to define for the 
highway traffic without the presence of the bottlenecks. 
It Is, however, unambiguous that when the traffic den¬ 
sity exceeds Pc, one almost never see wide moving jams 
evolve spontaneously from the free flow, unless something 
drastic (for example, an accident) happens. Instead the 
traffic breaks down to the congested (or the “synchro¬ 
nized”) phase which can last up to an hour; the wide 
moving jams, on the other hand, emerge from the con¬ 
gested phase via the “pinch effect” or merging of numer¬ 
ous narrow jams. This qualitative mechanism should be 
predicted by a useful model. 

The spatiotemporal patterns of the highway traffic 
with the presence of bottlenecks are also very Important 
In characterizing the traffic dynamics. The distinction 
between the “synchronized flow” and the wide moving 
jam can be unambiguously stated by looking at the be- 
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haviors of their respective downstream front. The down¬ 
stream front of the wide moving jam will move with the 
same characteristic velocity Vj when passing through the 
bottleneck, while the downtream front of the “synchro¬ 
nized phase” is pinned at the bottleneck]!]. When the 
bottleneck is in the form of an on-ramp through which 
additional vehicles are injected into the main highway 
traffic, a wide moving jam passing through the bottle¬ 
neck can either induce or suppress the congested flow 
upstream of the bottleneck, depending on the strength 
of the traffic flow along the main highway traffic Qm and 
from the on-ramp Qin ■ In general the congested flow can 
last for a few kilometers, and the wide moving jams tend 
to evolve from the congested traffic and move upstream. 

Detailed studies of the congested phase, including the 
“general pattern” and the “synchronized phase”, also 
show that when the bottleneck strength increases, the 
mean frequency of the moving jam emergence becomes 
greater, and the region of the congested traffic upstream 
of the bottleneck is also smaller. Another important ob¬ 
servation of the highway traffic indicates that there are 
still very significant fluctuations in the average velocity 
of the vehicles even in the region when the vehicle density 
is very high. These two fundamental empirical features 
were first pointed out to illustrate the inadequacy of the 
early GM model classes|l|, which seem to predict exactly 
the opposite behaviors. They are thus very important 
gauges in testing the usefulness of any constructed traffic 
models. The empirical features listed in this section is 
summarized in Table. |T| 


III. A GENERAL FRAMEWORK IN MODEL 
CONSTRUCTION 

A microscopic traffic model requires the understand¬ 
ing of the dynamics of a single vehicle based on the in¬ 
teraction with its environment. Unlike classical physical 
systems we cannot write down equations based on the 
symmetry of the system: there is no spatial or temporal 
symmetry in the traffic system, and even the interacting 
components are not identical to each other. In fact, the 
complex dynamics of the traffic system can be a result of 
the following (non-exhausive) factors: a), the non-linear 
interaction between individual vehicles; b). diversity of 
the vehicle/driver types; c). stochasticity of the driving 
behaviors; d). inhomogeneity of the traffic lanes; e). time 
dependence of the vehicle number and driving behaviors. 
To understand any empirical features observed in the real 
traffic, it is very important to show which one or few 
of those factors (and not others) are fundamentally re¬ 
sponsible, and this should be reflected in the constructed 
model. From a more theoretical point of view, one would 
also like to understand what interesting phenomena can 
result from the non-linear dynamics alone, independent 
of all other factors. 

In reality almost all highway traffics have multiple 
lanes, and the behavior of lane changing and overtaking 


can be important for certain empirical observations |321 
I33j . In this paper, however, we ignore multiple lanes 
completely by modeling the traffic as a one-dimensional 
system. Each component, or vehicle, is labeled by a sub¬ 
script n, which increases sequencially in the direction of 
the highway traffic, indicating no overtaking. We try to 
base our model on what we actually observe experimen¬ 
tally as much as possible. The acceleration of the 
vehicle is most generally given by 

O-n — ^ n,{si} ({fi}) (2) 

where {s^} and {ti} combined is the collection of envi¬ 
ronmental factors that influences the acceleration. The 
separation of these factors into groups {si} and {ti} is ar¬ 
bitrary, but from the modeling perspective {si} contains 
all the unimportant factors we would like to average over. 
This is because Eq.([^ by itself is not useful for analytic 
or numerical calculations. One can, however, repeatedly 
measure a„ over a wide range of {si} and {t^} and aver¬ 
age over {si}, which is formally represented as follows 

~ ^ ({^*}) ~ ({^*}) ( 3 ) 

where Nq is the proper normalization factor. If we also 
assume identical drivers in the traffic system, one also 
need to average over all the vehicles on the highway to 
obtain 

1 ^ _ 

k=l 

where N is the total number of vehicles. One can also 
give a time delay on the LHS of Eq. (lll-Eq.Q to model 
the reaction time of the drivers. Since the reaction time 
tends to be small and the acceleration is already the sec¬ 
ond time derivative, we will not consider it here. More 
importantly Eq. and Eq. can be empirically mea¬ 
sured after proper averaging. In practice, one does not 
need to know the details of {si} to obtain Eq.Q. A rea¬ 
sonably sufficient process is to record {a„, {ti}T of many 
vehicles for a long period of time over diverse environ¬ 
ments. For each set of {ti} one can average over the 
corresponding an to obtain Eq.Q. 

The choice of parameters in {ti} is motivated phys¬ 
ically and for convenience. Since we are interested in 
the spatiotemporal characteristics of the traffic dynam¬ 
ics, it is natural to take {ti} as a collections of posi¬ 
tions and velocities. The most intuitive choice is {U} = 
{hn, An„, Vn}, where is the bumper to bumper head¬ 
way of the vehicle, and At;„ = Vn+i — Vn is the 
velocity difference between the two consecutive vehicles. 
We are thus looking at the equation of motion of a one¬ 
dimensional system of identical components with a near¬ 
est neighbour anisotropic interaction. All other factors 
are averaged over and from Eq.Q we obtain: 

— fo Afyi, Vn) 


(5) 
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Flow-Density Diagram 

Emergence of the wide moving jams 

Congested traffic at 
an on-ramp bottleneck 

Pseudo-linear relationship 
when density is low 

The congested traffic (“synchronized phase”) 
can last up to an hour 

A wide moving jam passes 
through the bottleneck 
unaffected 

Large scattering of the congested 
flow-density data points 

Wide moving jams mostly emerge 
from the congested traffic 

A wide moving jam may induce or 
supprese congested traffic 
at the bottleneck 

The “hysteresis effect” 

“Pinch effect” and the merging of 
numerous narrow jams 

The region of congested traffic 
gets smaller with greater 
bottleneck strength 

Significant velocity fluctuation at 
very large vehicle density 


The frequency of the emergence 
of the moving jams increases with 
greater bottleneck strength 

Quantitative features: 

^max ? pdj 5 Pc ? Pj 5 Qc? Qdj 




TABLE I. The list of commonly observed and well-defined empirical observations from the German highway systems. 


where for notational convenience we remove the bar rep¬ 
resenting the average taken in Eq.([^ and Eq.Q. We will 
show in Sec. 0 that Eq.([^ contains the minimal set of 
parameters to capture the empirical features discussed in 

Sec. HU 

It is a data and labor intensive task to obtain empir¬ 
ically the exact model as defined in Eq.([^, but we do 
know from common driving experience that a„ should 
increase with increasing , Ar;„, but decrease when Vn 
increases. In addition, we make two key assumptions. 
We first assume there exist solutions to /o {h, 0, v) = 0; 
the solution is basically a statement that for = ho, 
Avn = 0, all vehicles are equally spaced apart travel¬ 
ing at the same velocity with no acceleration. We de¬ 
fine this as the ground state of the traffic system at av¬ 
erage density . This gives the implicit solution(s) 

/o (jio,0,Vop^ (^o)J = 0. For each hg there can exist 

f A;') 

more than one Vop ^ indexed by k, leading to more than 
one ground states with different velocities at the same 
traffic density. 

The second assumption is that /o is smooth around 
Avn = 0 and Vn = Vop'^. The two physical scales of 
the traffic system is pj, the maximum vehicle density 
which occurs within a wide moving jam; and Vmax, the 
maximum velocity. Using the dimensionless quantities 

kji = hjiPj ,Vn — '^nt^maxi ^op — ^op j^raax and = 

AvnlVniax, we define 

/o = Pj^K^fo (jln, AVn, (6) 

where fg is also dimensionless and k = pjVmax- The 
assumption allows us to do Taylor expansion around each 


ground state as follows 



AVn—0 
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+ 0(2) 


(7) 


when Vn — Uo^p ^ (jin^ as well as Ai;„ are small, and where 
0(2) contains terms of higher orders of expansion. Here 
we also define a„ = anPj, so the only dimenional scale in 
the equation is k which gives the inverse time. 

A few comments are in order here. The averaging 
process performed in Eq.([^ and Eq.Q leads to a time- 
independent, deterministic model with identical drivers. 
These are the general assumptions for most microscopic 
traffic models in the literature. The resulting model in 
Eq.§ can be easily generalized to more (long-ranged) 
interactions, for example by including the next nearest 
neighbour {vn+2 G {U})) or backward looking {vn -i € 
{U}), though we will show in Sec. 0 and Sec. they 
are not necessary. 

On the other hand, it is unlikely that Eq.([^ will di¬ 
verge in the limit /i„ —>■ 0. In fact since Eq. (§ is in 
principle averaged from the empirical data, a„ must be 
bounded both from the above and from below. However, 
popular models like the IDM with an artificial divergence 
when the headway goes to zero also can be expanded 
around its unique ground state at each average density. 
One would in fact expect any model with a fundamen¬ 
tal diagram to have an equivalent model in the form of 
Eq.§, with a unique optimal velocity function that can 
capture the same physics of the traffic dynamics. 
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The “speed adaptation model” in proposes two dif¬ 
ferent optimal velocity functions in two velocity ranges. 
If one assumes /o is analytic, this corresponds to expand- 
ing around two different Vop . One should note that 
if we assume smoothness of fo and if there exists mul¬ 
tiple Vop\ we would expect dfo/dvn > 0 when evalu¬ 
ated at Avn = 0,Vn = Vop^ for some k. Around these 
points Eq.Q is generally unstable with both velocities 
and headways diverge over time. Consequently the low¬ 
est order approximation in Eq.Q becomes invalid and 
most probably the system will settle into one of the sta¬ 
ble regions where dfo/dvn < 0. Thus if experimentally 
fo is found to be analytic with multiple Vfp for certain 
range of /i„, then the “speed adaptation model” can be 
justified microscopically. 

Most of the complicated “three-phase” microscopic 
models proposed assumes that fo (hn, Avn, Vn) is not an¬ 
alytic when hn is smaller than the so-called “synchro¬ 
nization gap”. Unlike the IDM models where the di¬ 
vergence of the acceleration is purely an artificial mod¬ 
eling tool, the non-analyticity of fo can in principle be 
checked with the experimental data. It would be interest¬ 
ing to see if the exact fo from the empirical measurement 
shows non-analytic behavior. In this framework the non- 
analyticity of fo is the necessary condition for the mul¬ 
titude of steady states with a non-unique relationship 
between the flow and density. 


IV. EMERGENT PROPERTIES OF THE 
ANALYTIC MODEL 

We now focus on Eq.Q with one particular Vop^ 
and discuss the mathematical properties of the result¬ 
ing coupled ODEs with non-linear interactions. We first 
rewrite Eq.Q in a simpler form, ignoring higher orders 

of (h„ - Vof) : 

dn = Ki (jln^ (V}p'> (jln^ - Vn^ + g (jln, AVn'j ( 8 ) 

here ki (Jin^ = K^dfo/dvm and the second term on the 
right still keeps all the higher orders of Ah„. This is 
just a general form of the optimal velocity (OV) model, 
and its mathematical properties are quite well known. 
In addition to briefly discussing those properties, we will 
also introduce an emergent intrinsic length scale, which is 
not only theoretically interesting by itself, but also serves 
as an important tuning parameter for Sec. |Vj 

To simplify the discussion we use the simplest form 
of the OV model as an example, and all the relevant 
properties are universal and can be qualitatively applied 
to Eq.([^ unless otherwise stated. In this simplest case 
we ignore the dependence of the acceleration on Au„ and 
choose the most popular optimal velocity function[5]: 

UjO = Ui + U 2 tanh (Ci- C 2 ) (9) 


We also take the special case ki 


(hn'j 


= Xk^, indepen¬ 


dent of hn- Dehning = Cihn — C 2 , and rescaling the 
time t —>■ KCiV 2 t, a simple transformation gives us the 
equivalent form 


Sn + KoSn = Ko (tanhs„_|_i — tanhs„) 


( 10 ) 


with the only dimensionless parameter ko = X/ (C'iV 2 )- 


We will now focus on Eq.(10l, where s„ is dimension¬ 
less. The change of variable and the scaling away of 
the dimensions not only tells us that seemingly different 
driving behaviors are actually equivalent, it also makes 
the symmetry of ODE’s in Eq.p^ explicit. One should 
note by definition a physical which is always positive, 
can lead to negative s„ depending on the parameters in 
Eq. (|^. Linear analysis leads to a stable phase of s„ = sq 
against small perturbation, and the spinodal line (or the 
neutral stability line) is given by 


2secn So = ^o- 


( 11 ) 


In the regime |so| > Sd = |sech“^ a/k/ 2|, a small pertur¬ 
bation to a uniform headway sq with s„(t —>■ 0) = so-l-(5s„ 
leads to s„(t —)■ 00 ) = So, so this regime is linearly sta¬ 
ble. The uniform headway solution is the ground state 
defined in Sec.|III[ Note Eq.(ll I is only exact in the limit 


when the perturbation goes to zero; close to the spinodal 
line, the uniform headway configuration is metastable, a 
large enough perturbation will also lead to the formation 
of clusters. 

We now show that the coexistence curve that separates 
the metastable phase and the absolutely stable phase 
can be numerically obtained from the cluster structure. 
Eirstly, in the regime |so| < Sci, it is well known that 
small perturbations will grow in time with the formation 
of clusters, as shown in Eig.([^, where a random initial 
condition settles into a conhguration with the majority 
number of vehicles having two extremum headways given 
by ±Sc 2 . As smaller s„ implies higher physical vehicle 
density, vehicles with headway — Sc 2 form clusters or jams 
of very high density with minimal velocity, while vehicles 
with headway Sc 2 moves with very high velocity, forming 
anti-clusters. Interestingly like Sd, the numerical value 
of Sc 2 only depends on kq but not on sq) so the cluster 
structure is unique once the parameters in the model is 
fixed. 

Secondly the number of vehicles involved in the “kink” 
or “anti-kinks” are independent of sq and the total num¬ 
ber of vehicles N. A “kink” is the “go front”, or the 
transition region from a cluster with ~ — Sc 2 to an 
anti-cluster with s„ ~ Sc 2 , while an ” anti-kink” is the 
“stop front”, or the transition region from an anti-cluster 
to a cluster. Thus for large N we can ignore vehicles in 
the “(anti-)kink”, and the number of vehicles in the clus¬ 
ter is given by 


^ 2 


N Sc2 - So 


Sc2 


( 12 ) 


Clearly for sq ^ Sc 2 : no clusters can be formed, given 
random initial perturbations of any magnitude, as long 





7 


Mathematical Expression 

Assumptions Implemented 

Comments 

dfi — ^ ri,{si} (^717 Vn-, ^Vri) 

ho , Vo and Avo chosen 
as the important parameters 

More (or different) parameters can be 
chosen as important, leading to different 
empirical features captured or lost by 
averaging. While T is formal, the exact 
form of fo can be measured experimentally 

dn ~ fn {,^ni '^717 ^'^71) 

The non-essential parameters in {si} 
can be averaged over 

d-n — fo (^77,7 '^nt ^'^n) 

Identical Drivers 

More than one species of 
vehicles can be included 

/o(/io,0,vip' {ho)) =0 

At least one ground state exists 

Each of the assumptions 
can be experimentally verified 

(^") (hn)) 

fo is analytic around those solutions 

dfi = (^77) (^77 

Unique ground state 
at each average density 
(2-phase models with 
a fundamental diagram) 

a„ = no {h„) {vn - Vop {h„)) + g (/i„, Aii„) 

Keeping the lowest order 
of expansion around Vop 

Higer orders can be ignored 
if Vo — Vop is small 

compared to Vmax in the congested phase 

a„ = Ko {vn - Vop {ho)) + g {h„, Av„) 

Coefficient kq independent of ho 

This assumption can be 
experimentally verified 

g {ho, Avo) = XiAvo + A 2 Av„ 

A particular form of g is chosen 

This is a mathematically convenient 
form to include the non-linearlity of g {Av„) 


TABLE II. Various stages of assumptions implemented for the construction of the effective model in Sec. with their corre¬ 
sponding mathematical expressions. 


as = 0- Similarly, no anti-clusters can exist for 

So < — Sc 2 - We thus identify Sc 2 as the coexistence curve 
and plot it together with Sd in Fig.Q. The 
numerically calculated coexistence curve and the spin- 
odal line coincides at the critical neutral stability point 
located at sq = 0, /c = 2, agreeing with the previous 
analysis [T7115^ . Note that s„ can be negative, and the 
physical vehicle density is calculated from the model pa¬ 
rameters Cl and C 2 . There is also a duality between 
So O —So, where clusters at so corresponds to anti¬ 
clusters at —So, and all behaviors at So are identical to 
those at — sq. This symmetry is entirely due to the fact 
that the RHS of Eq.(10l is odd. In the more general 
model of Eq.([^, this duality can be broken and it is thus 
not universal. 

Progresses have been made in treating non-linear ODE 
describing car-following models analytically [THl [S21 EHl 
l60] : For Eq.(lO) it is generally accepted that one can 


do a controlled expansion near the critical neutral sta¬ 
bility point and close to the neutral stability line; the 
former leads to the modified KdV equations plus cor¬ 
rection terms, that gives the approximate “(anti-)kink” 
solutions; the latter reduces the original model to the 
KdV equations plus corrections that give rise to soliton 


solutions|57j. However, away from the neutral stability 
line, it is clear from the numerical calculation that if one 
makes the vehicle index continuous, the transition be¬ 
tween the two extremum headways is discontinuous and 
analytically intractable. 

One can, however, show that the “kink” and “anti¬ 
kink” of a single cluster move at the same velocity, by 


taking s = Yli=i 


For the “kink”, the i**' vehicle is 


located in the cluster, while the vehicle is located in 
the anti-cluster. From Eq.(lO) we have 


s + KqS = 2kq tanh 3^2 


(13) 


The relevant set of solutions is s = 2tanh(sc2)i + C, 
where C is an unimportant constant of integration. This 
gives the velocity of the “kink” as the number of vehicles 
per unit time as follows 


Vk = 


tanh Sc 2 
Sc2 


(14) 


The velocity of the “anti-kink” is calculated similarly, 
thus Vk gives the velocity of the cluster, which again is 
independent of the vehicle density of the traffic lane. Here 
we make the assumption that for vehicles far away from 
the “kink” or the “anti-kink”, their headway takes the 




















Car Index 



So 


FIG. 2. The plot of the headway as the function of the vehicle 
index, when a jam or a cluster is formed. This cluster configu¬ 
ration evolves from a random initial headway distribution, as 
shown in the top inset. The bottom inset is the spinodal curve 
(the solid line without circles, plotted from Eq.( 111), and the 
coexistence curve from the numerical calculations (the solid 
line fitting the solid circles). The solid circles are numerically 
observed extremum headways at different kq. 


FIG. 3. (Golor online) The probability of having one to six 
clusters in a single traffic lane, plotted as the function of the 
initial average headway sq. The probability is calculated with 
three hundred vehicles and random initial headway perturba¬ 
tion. Inset: The probability of having only one cluster, as 
the function of the initial average headway. The probability 
is calculated for three hundred vehicles to seven hundred ve¬ 
hicles, showing that in the limit of large number of vehicles, 
the probability curve converges to a well-defined limit. The 
probability is calculated at t = 30000s. 


value of ±Sc 2 . More importantly, if we concatenate two 
clusters together, as long as the assumption holds (e.g. 
when the two clusters are far away), they will move at 
the same velocity and will never merge. 

One interesting universal aspect of the OV models is 
the non-trivial probability distribution of the number of 
cluster formations, when the ground state is randomly 
perturbed. One might expect that a random initial state 
like the inset of Fig.([^ should lead to a random number 
of clusters |?T]. at least in the limit of large N, subject¬ 
ing to the constraint of Eq.( 12). However, our numerical 
results show that the probability distribution of the num¬ 
ber of clusters is not random; it strongly depends on the 
initial headway Sq and Kq. An example of the proba¬ 
bility distribution is calculated and presented in Fig.(|^, 
by fixing the strength of the initial random perturbation 
and kq in Eq.(lO), and only vary the initial headway sq. 
For each value of Sqj sufficiently large number of random 
initial states are generated until the probability for each 
number of clusters converges. 

A few comments are in order here before we give a 
proper explanation. In Fig.(|^ we only plot the part 
where Sq is negative, because the probability distribution 
is identical for sq and —sq. For |so| > 0.87 we can see the 
final state is dominated by the one-cluster configuration, 
and this is true even for an infinitely long traffic lane 
as > oo; in this case, most probably one very large 
cluster develops, instead of several clusters with smaller 


lengths. As |so| decreases, the probability of having more 
than one cluster increases, and for |so| < 0.82, it is al¬ 
most impossible to have just one cluster. As |so| further 
decreases towards zero, the average number of clusters 
most probably will tend to infinity. This cannot be ob¬ 
served numerically for a finite number N, since at Sg = 0 
the total number of vehicles in the clusters is ~ A^/2(see 
Eq.@). 

To understand the probability distribution of the 
number of clusters, we characterize quantitatively the 
strength of interaction between two clusters by the time 
it takes for them to merge. It is useful to plot dsn/dt in¬ 
stead of Sn as a function of the vehicle index n. The 
“kinks” and “anti-kinks” lead to exponentially local¬ 
ized “quasisolitons” of opposite charges (see Fig.([^), 
which closely resemble the “autosolitons” in dissipative 
non-linear svstems[54j. When quasisolitons of oppo¬ 
site charges annihilate each other, two clusters or anti¬ 
clusters merge into one. We numerically observe that the 
time needed for annihilation, ta, increases exponentially 
with the number vehicles n between the peaks of these 
two quasisolitons, giving the relationship 

ta - (15) 

One thus note that when |so| increases, the cluster (for 
So > 0) or the anti-cluster (for sq < 0) region gets nar- 
rower(see Eq. iHl)), leading to higher probability of short 
distances between quasisolitons. Thus the probability of 
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FIG. 4. The quasisoliton structure of ds„/dt as a function 
of the vehicle index. The fewer the vehicles involved in the 
quasisoliton, the smaller the width of the quasisoliton, which 
depends only on kq and not on the initial headway sq. By 
convention a kink gives a positively charged quasisoliton as 
shown in this figure. An anti-kink gives a negatively charged 
quasisoliton. The top inset shows the dependence of the an¬ 
nihilation time ta on the number of vehicles between the qua- 
sisolitons of opposite charges, the exponential fit is numeri¬ 
cally perfect. The intrinsic scale as a function of Sc 2 is shown 
in the bottom inset. 


<o 


m 


® CM 
< 



0.74 0.76 0.78 0.80 0.82 0.84 0.86 

So 


0) ^ 
oi 


< CM 


1.50 


1.55 


1.60 


1.65 


1.70 


Sc2 


FIG. 5. The average number of clusters of a single lane traffic 
as a function of the intial headway (the top plot, while keeping 
the perturbation strength and kq fixed), and as a function of 
Sc 2 (the bottom plot, while keeping the perturbation strength 
and So fixed). 


having multiple (anti-) clusters is suppressed, as shown 
in Fig.([^. The intrinsic “scale” uq in Eq.( 151 depends on 
Sc 2 or Kq, which is also plotted in Fig. This is analo¬ 
gous to the interaction and collapsing of kinks and anti¬ 
kinks in the Ginzburg-Landau theory [KH]. though here 
the total number of vehicles in the cluster has to sat¬ 


isfy Eq.(12|, so that at least one cluster will remain for 


a finite system with periodic boundary condition. Thus 
the greater the intrinsic scale, the stronger the interac¬ 
tions between the quasisolitons, so this scale can be used 
to quantify the absolute value of the quasisoliton charge. 
The interaction leads to merging of clusters, reducing the 
probability of having multiple clusters in the traffic lane. 
Fig.@ will look qualitatively the same if the x-axis is 
replaced with increasing Sc 2 - The dependence of average 
number of clusters as a function of sq and Sc 2 are plotted 
separately in Fig.([^, numerically supporting the above 
explanation [53j . 

One should note that the symmetry Smax = —Smin = 
Sc 2 is a result from the fact that the RHS of Eq.(lO) 
is odd in s„. This symmetry forbids us to tune the 
corresponding h^ax and hmin independently. The sym- 
metry could be broken when Vop ^, or additional terms 
from g (hn, Az;„) in Eq.(|^ contains components that are 
even in s„. In these cases, all three emergent quantities 
Smax.Smin (or the physical headways hmax,hmin) and 


Uq can be tuned freely with the parameters in the model, 
providing the necessary degrees of freedom in Sec. |V] to 
capture the empirical features in Sec. |ll] 

Understanding the multicluster solutions and the phys¬ 
ical significance of the intrinsic scale Uq is important in 
explaining some of the essential features of the traffic dy¬ 
namics. While spatially random perturbations grow and 
interact with each other in the metastable and linearly 
unstable region and eventually form very wide clusters 
corresponding to the large moving jams, physically no 
characterizes the time scale over which the intermedi¬ 
ate transient states can last, as well as the width of the 
narrow jams that can be detected by the actual measure¬ 
ment. Numerical calculation also shows no can be used to 
tune the maximum acceleration that would normally oc¬ 
cur in the congested traffic (which is much smaller than 
the physical limit of the acceleration of the actual ve¬ 
hicle), making it an important parameter in tuning the 
effective model. 


V. ALGORITHM FOR TUNING THE TRAFFIC 
MODEL 

We will now proceed to construct the simplest effective 
model, or the minimal model for the real traffic dynam¬ 
ics, that can capture what we observe in Sec. The dis- 
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cussions in Sec. |IV| is universal, and for a general model 
given by Eq.([^, there are also three emergent quanti¬ 
ties hjnaxjhmin ^nd uq. Physically, Uq also controls the 
maximum acceleration in the congested traffic, and the 
interaction between clusters is crucial for the time scale 
of the evolution of the wide moving jams. We make the 
following assumptions to start with the simplest possible 
model from Eq.Q: 

a) Only the linear order of (vjp^ — Vnj is kept. The 
model is thus reduced to Eq.([^. 

b) In Eq.(|^ we make both ki and g independent of 
the headway /i„, so the only headway dependence 
is within Vop . 

c) We assume /o decreases monotonically with respect 
to increasing for any at At;„ = 0, thus there 

(k) 

is only one Vop . 

Following assumption c) we remove the subscript (fc) in 
Vop. Assumptions in a) and b) are simplifications of the 
exact model, which can only be justified if — Vop is 
small compared to Vmax, and both ki and g depend 
weakly on throughout the time evolution of the traffic 
dynamics. If we only keep the linear order in Av^ in our 
expansion and the expansion coefficient independent of 
hn, the resulting model is the full force velocity model[6]. 
It does not however contain enough degrees of freedom 
to capture all the empirical features in Sec. |TT] because 
of the symmetry Smax = —Smin (see Sec. |IV[ ). One can ei¬ 
ther make the coefficient of expansion dependent on /i„, 
or keep the higher orders in For simplicity we choose 
the latter option. Here we postulate the empirical fea¬ 
tures in Sec. can be universally captured by the three 
emergent quantities of the general traffic models; the sim¬ 
plifications we undertake only remove non-essential mi¬ 
croscopic details we are not prepared to capture. 

We use the well-studied A-5 North German High¬ 
way from Kerner|2T] as an example. Tuning the model 
only requires the information from the flow-density plot, 
which consists of the free flow part (where the flow de¬ 
pends approximately linearly on the density), the con¬ 
gested part (with a collection of randomly scattering data 
at higher density with suppresed flow) and the wide mov¬ 
ing jam given by the “J line”. The list of statistically 
robust quantities from the flow-density plot we use are: 

Emax = lim dQ/dp ~ 42tos“^, Qmax ~ 3000veh/h, 

p—^0 

Pc ~ 30veh/km, Qdj 2000veh/h, 

Pdj n.bveh/km, pj ~ 125veh/km (16) 

Here pc is the critical density of the highway at which 
Qmax) the maximum flow, is observed. Qdj and pdj are 
the flow and density downstream of the wide moving jam 
respectively, while pj is the density within the jam. For 
the lack of the raw traffic data, all the numerical values 
are rough estimates only, and for our purpose of illustra¬ 
tion that is sufficient, as we do not need to fine-tune the 


model to simulate the qualitative empirical features. We 
also assume on average the length of the vehicle Q = 5m, 
and by identifying the parameters of the cluster structure 
with the characteristic parameters of a wide moving jam 
we have the following relationship: 


Vop (oo) — Vmax^ Vop {hmax^ — '^dj — Qdj / Pdj 

bop {hmin') — O5 kjjp (her) — '^cr — I?max/Pc (^^) 


The two other characteristic velocities from the flow- 
density plot are Vj = Qmax/ {Pc ~ Pj)i the velocity 
of the downstream front of a wide moving jam, and 
V(j; — (Qmax Qdj} / (per Pdj}-) the velocity of the 
downstream front between Qmax and Qdj. The clus¬ 
ter parameters are given by hmax = Pd^ - 4, = 


P 3 


4: her — Pc 4- 

The simplest form of Vop has been suggested in[i 


for its analytic tractability; it however does not capture 
the correct funda men tal diagram in the free flow phase. 
Here we solve Eq. ( |17[ ) most simply with a piecewise func¬ 
tion passing throu^(hmm, 0), (hm, Ver}) (hmax, Vdj} and 
bounded at Vmax) so as to fix the quantitative features 
of the real traffic dynamics[5^ in the free flow phase. 

- hmax} - Ic we have: 


Eop (h) = 


Vc; 

r-Vdj 

0 



Ver 

hcj 



zr-Vdj 

her 

• — k-max 


Em 


h ^ h-foin 

{h + Zc) T E,' h cr > h> h 

min 

{h -1- Zc) -t- Vc he ^ h ^ her 
h > he 


(18) 


Thus the fundamental diagram is defined by Eop (see 
Fig.@), capturing the quantitative features of the em¬ 
pirical flow-density plot. The next step is to tune kq 
and g {Avn} so that the cluster solutions have the de¬ 
sirable hmax)hmin and tiq. We adopt the reasonable as¬ 
sumption that the maximum acceleration for the vehi¬ 
cles in the stop-and-go wave should be within the range 
of ±3ms“^. Small variations around this quantitative 
assumptions do not qualitative alter the arguments and 
conclusions in this paper. Given that three parameters 
need to be fixed, g {Avn} has to be non-linear and con¬ 
tain terms that are even in Aw„ to break the symmetry of 
Smax = —Smin (see Sec. IV). We choose to simply adopt 


the AFVD model[S] (in the case where A 2 ^ 0) as follows 


9 {Avn) = XiAvn -H A2|A?;„ 


(19) 


The effective model is now completely defined, and with 
numerical calculations the fitted parameters are kq = 
0.1s“^, Ai = 6.2 and A 2 = —2.9, corresponding to hmin 
3m and hmax ~ 52m. 

A summary of our algorithm is in order. The opti¬ 
mal velocity function we chose defines the fundamental 
diagram: it gives the correct maximum average veloc¬ 
ity of the traffic system (when the density of the traffic 
on the highway is very small). It also gives the right 
{Pdj)Qdj} and (pcjQmax) paus on the flow-density plot, 
where the traffic is still in the free flow phase. In ad¬ 
dition, it gives the maximum density pj of the traffic. 
While these are just some of the special points on the flow 
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density plot, the tuning of the other parameters in the 
model (ko,Ai,A 2 ) makes sure {Qdj,Pdj) corresponds to 
the characteristic flow and density downstream of a wide 
moving jam, and pj corresponds to the density within 
a wide moving jam. In addition, it also makes sure the 
acceleration of vehicles in congested traffic is physically 
reasonable. 

The parameters in the model all have very clear physi¬ 
cal meanings; the general features of the optimal velocity 
function is also quite intuitive. Both Vop (hn) and g (Ai;„) 
are not smooth, which is not physical. They are, how¬ 
ever, just unimportant artifacts of the model that can be 
easily (but tediously) removed mathematically without 
affecting any of the conclusions or the predictive pow¬ 
ers of the model. One should note that the optimal ve¬ 
locity function and the parameters in the model do not 
explicity tell us about the transition from the free flow 
to the congested flow, or the maximum flow that can be 
achieved, as well as the qualitative features of the com¬ 
plex spatiotemporal patterns. All these features will be 
predicted by the dynamics of the model which we will 
show in the following section. 


VI. PREDICTIONS OF THE EFFECTIVE 
MODEL 


We now proceed to examine what the minimal model 
predicts about the traffic dynamics. Given the param¬ 
eters in the model, the free flow in the stable phase 


is given in the region p < {hmax + 0 


-1 


llveh/h. 


The metastable region is given by 17veh/km ^ p 
30veh/km. In this region, a large enough perturba¬ 
tion will grow in time and leads to the instability of the 
free flow and formation of the jams. The empirical fea¬ 
ture that the free flow persists up to the critical density 
Pc ~ 30veh/h is nicely predicted by the fact that for the 
metastable region with density smaller than ^ 30veh/h, 
the perturbation needs to be greater than the average 
vehicle headway for the free flow to be unstable (see 
Fig.®), which is unlikely without collisions. Thus the 
effective model captures Qmax and pc quite accurately, 
even though Eq.(18) in no way guarantee the stability 
condition agreeing with the empirical data. 

To study the congested traffic and the evolution of the 
wide moving jams, both periodic boundary condition of a 
single homogeneous lane and open boundary condition of 
a single lane with the presence of an on-ramp bottleneck 
are simulated. For periodic boundary conditions, the ini¬ 
tial condition is chosen to have random fluctuations of 
the headways with different average density. Though in 
the long time limit wide moving jams eventually form 
for average density p > 29veh/km, the intermediate pro¬ 
cess can be quite complicated. When the average density 
is very close to the phase boundary, or the coexistence 
curve given by {hmax + 0 ; very large perturbations are 

needed to nucleate a wide moving jam via the well-known 
“boomerang” behavior |H when the average density 



FIG. 6. (Color online) Characteristics of the effective OV 
model. Left: The piecewise optimal velocity function. Left 
inset: The magnitude of perturbation {Ah) needed to form 
cluster solutions in the metastable region. Cases for the sin¬ 
gle vehicle perturbation (red) and the random perturbation 
(black) are plotted. Right: The actual flow-density diagram 
including both the free flow and the congested flow. The plot 
is obtained from the numerical calcuation with open bound¬ 
ary condition and an on-ramp bottleneck. The left vertical 
line gives pdj and the right vertical line gives per- 


increases further, the “pinch effect” [3] is observed at mul¬ 
tiple locations leading to multiple narrow jams; at high 
density numerous narrow jams form relatively quickly, 
and over time these narrow jams interact and merge into 
a few wide moving jams (see Fig. 0). 

While the “pinch effect” and the formation of numer¬ 
ous narrow jams are well-known in the literature [3], the 
existence of the “boomerang” behavior is still debated [31 
[2Hj . The absence of “boomerang” behavior may also 
due to the rarity of very large perturbations on a multi¬ 
lane highway, when accidents and bottlenecks are absent. 
Even for moderately large perturbations, it takes more 
than an hour for a closed traffic system to develop the 
“boomerang behavior”. One should also note that based 
on the numerical simulation, it takes 30 ^ 60 minutes for 
the wide moving jams to eventually emerge from a ran¬ 
dom initial condition via complex intermediate states. 
Thus comparing numerical results with a fixed number 
of vehicles and periodic boundary condition to the empir¬ 
ical observations can be extremely tricky. The real world 
traffic, being an open system, does not maintain its ve¬ 
hicle density and the total number of vehicles over an 
extended period of time; variation of the average density 
within the metastable/unstable region leads to a mix¬ 
ture of long lasting intermediate states, numerous narrow 
jams and occasional wide moving jams. 

From both the empirical validation and transportation 
engineering points of view, numerical simulations of the 
effective model in the presence of bottlenecks are more 
crucial. In our simulation with open boundary condition 
the virtual sensor measures the flow and average velocity 
of the passing vehicles in exactly the same way as the traf¬ 
fic sensors installed in the real world highways [27] . We 
use the idealized initial condition with a constant main 
traffic flow Qm and an on-ramp flow Qi-n- At low enough 
Qin the free flow is maintained, though in the linearly 
unstable region the on-ramp flow has to be close to zero. 
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When Qin increases, corresponding to the increase in the 
bottleneck strength, the congested flow develops immedi¬ 
ately upstream of the bottleneck. This region with length 
Lc share the characteristics of the “synchronized flow” 
(see Fig. ®,Fig.0). Narrow jams form upstream of the 
congested flow, and wide moving jams appear upstream 
of these narrow jams from merging of the narrow jams 
and the “pinch effect”. The “boomerang” behavior is 
also observed (see Fig.([^)). The congested flow can be 
either spontaneous or induced by a passing wide moving 
jam. In general, L^. can be as long as Akm and decreases 
with the increase of Qin- 

Previous GM based models were criticized in [i] based 
on some fundamental empirical observations of the con¬ 
gested phase at the bottleneck, as well as on the absence 
of the homogeneous congested traffic (HOT) empirically. 
In contrast, the effective model we constructed agrees 
with the empirical observation that increasing the bot¬ 
tleneck strength leads to higher frequency of moving jam 
emergence and smaller Lc (see Fig.Q). In fact this is the 
most common situation for various different Qm- Numer¬ 
ical calculations also show the metastable phase in the 
high density region is very narrow. The model actually 
predicts complicated spatio-temporal structures for the 
congested traffic at very high density, with traffic flow 
fluctuate between zero to 500 veh/h, agreeing with the 
empirical observation in[4]. This is simply because small 
perturbation is linearly unstable even in the region of 
vehicle density up to ~ lOOveh/km. 


VII. SUMMARY AND OUTLOOK 

In summary, we have presented a systematic way of 
constructing traffic models that in principle can capture 
all the empirical features for which stochasticity as well 
as diversity in drivers and vehicle types do not play a 
fundamental role; only the averaged non-linear interac¬ 
tion between the vehicles is fundamental. The exact form 
of the traffic model given by Eq. <§ can be obtained by 
the empirical measurement. In the case that the model 
is analytic, one should obtain various effective models by 
expansions around the ground states defined in Sec. m 
The compromise between keeping the model simple and 
capturing more microscopic details can now be done in 
a systematic way by choosing the appropriate set of pa¬ 
rameters to average over, and by making various approx¬ 
imations in the expansions. 

In addition, we proposed a simple algorithm to jus¬ 
tify the approximations we made by expanding the traffic 
model around a unique ground state. The resulting min¬ 
imal effective model shows that the physics of many em¬ 
pirical observations of the highway traffic dynamics can 
be captured by a deterministic effective model based on 
a simple optimal velocity function we proposed. In this 
framework the congested traffic is characterized by long 
lasting transient states of the model, from which the wide 
moving jams evolve from the “pinch effect” or the merg- 



FIG. 7. (Color online) Various traffic patterns from the min¬ 
imal effective model with the velocity color plot. Plot a)~c) 
are numerical calculations from the periodic boundary con¬ 
dition; the rest are from the open boundary condition with 
the presence of an on-ramp bottleneck, a). The “boomerang” 
behavior with average density close to the coexistence curve, 
which starts at t =0 min at around distance ~ 16 km. b). 
The “pinch effect” in the metastable region, c). Numerous 
narrow moving jams in the linearly unstable region, d). A 
wide moving jam followed by the congested traffic at the bot¬ 
tleneck when the main traffic density is metastable, e). The 
formation of the congested traffic at the bottleneck when the 
main traffic density is stable, f). Small Qin with metastable 
main traffic flow, a wide moving jam effectively stops the con¬ 
gested traffic because the density downstream of the jam is 
lower than the that of the main traffic, g). A wide moving 
jam passes through the bottleneck with its downstream front 
maintaining its characteristic velocity, h). A wide moving 
jam induces congested traffic at the bottleneck. The main 
traffic density is lower than the density downstream of the 
wide moving jam. 


ing of the narrow jams. Interestingly, our results imply 
it is probably difficult to distinguish between real traffic 
dynamics and those from identical autonomous vehicles 
with simple driving rules based on the empirical data, 
including the flow-density plot and the congestion pat¬ 
terns near the bottlenecks. One should note that even de- 
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FIG. 8. (Color online) Additional features predicted by the 
minimal effective model in Sec. [V] with the velocity color plot 
except the bottom right one. Top left: Congested pattern at 
the bottleneck with Qm = 18Q0veh/h and Qin = 360veh/h. 
Top right: Congested pattern at the bottleneck with Qm = 
1800veh/h and Qin = 50Qveh/h. Both the narrow jams and 
the wide moving jams emerge closer to the bottleneck, and 
fluctuating with higher frequencies. Bottom left: Dynamics 
of a model state with the wide moving jam adjacent to the ho¬ 
mogeneous congested state. The velocity of the downstream 
front of the wide moving jam is unaffected by the traffic con¬ 
ditions downstream of the jam. Bottom right: The flow time 
series when the traffic density is 80veh/km 

terministic models can simulate seemingly unpredictable 
dynamics because in reality the initial condition can be 
random and the highway traffic is an open system. 


While the minimal effective model we constructed give 
convincing evidence that many complex traffic dynamics 
can still be explained solely by the simple interactions 
between nearest neighbour vehicles, we would not ex¬ 
pect the model to explain all the interesting empirical 
features [3]. With more advanced techniques in collect¬ 
ing and analyzing empirical data, more empirical fea¬ 
tures can be properly defined, and we would expect to 
relax certain assumptions in our effective model(see Ta¬ 
ble. 0 to capture those features. One should note that 
the optimal velocity model class emerge naturally from 
an analytic traffic model, though such model can still be 
complicated if there is more than one ground state for a 
fixed average density. This, together with the possibility 
that the model could be non-analytic within certain den¬ 
sity range, can justify the parameter rich “three-phase 
models” from a more fundamental ground. A long term, 
microscopic measurement of the vehicle acceleration as 
a result of the interaction between its close neighbours 
under diverse environmental settings is currently work in 
progress. The empirical data, together with proper aver¬ 
aging, should be able to give us Eq. and thus provide 
insight on which effective model is more consistent with 
the nature. 
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